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ABSTRACT 

The commonly-used "4/3-Earth Nodel", in determining 
radar altitudes, has recently been shown to over-estimate 
curvature at elevations greater en several kilometers. 
It is further limited in that provisions are not made to 
account for commonly-occurring super-refractive atmos- 
pheric layers. Based on these shortcomings, a new model 
is constructed which incorporates an exponential refractiv= 
ity model proposed by recent investigators, together with 
superimposed low-level trapping layers of various inteén- 
sities. From this model, tables of heizhts are computed 
for altitude-determining radars. These heishts should be 
considerably more representative than those based on the 
'4/3-Earth Model. 

The author wishes to express his appreciation for the 
assistance oe Professor F. L. Martin and Mr. E. N. Ward 
of the U. S. Naval Postgraduate School in this investiga- 
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SYMBOLS AND ABBREVIATIONS 


Radius of the earth 

Effective radius of the earth 

Angle of ray emission, relative to the horizontal 
Ray-entrance angle, relative to the horizontal 
Ray-emergence angle, relative to the horizontal 
feight of trapping layer 

Thickness of trapping layer 

Refractive index of atmosphere 

Refractivity of atmosphere 

Height, measured from surface of the earth 
Range 


Relative curvature of ray 
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1. Introduction 

Atmospheric refraction of radar rays can seriously 
affect the accuracy of surface-based heisht-finding radars: 
The purpose of the paper is to consider some models of 
anomolous atmospheric refraction, and determine corrections 
necessary to observed radar altitudes in order to provide 


an accurate target altitude. 


Figure 1. Typical radar ray path. 
Consider.a transmitter, -T,.in.Figure »l,.,and .its_target, S. - 
Because of atmospheric refraction, the radar ray path will 
be, in general, a ,curved arc,.TS...iIn.order to determine 
accurately the altitude of the target, S, we must be able 
LO predict correctly the.curvature of the are TS along«all 
points.of the arc. 

In 1933, Schelling, Burrows, and Ferrell L4], set 
forth the effective earth's radius model to account for 
the effects of atmospheric refraction, better known as the 
"4/3-Earth ilodel", This is to say that if we were to depict 
the earth as a sphere of some radius suitably larger than 
that of the real earth, then normally downward bending 


radar waves could be drawn as straight lines. The accepted 
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standard effective radius has become 4/3 the radius of the 
earth. 

Recent investigations by Bauer, Mason, and Wilson eae 
of Lincoln Laboratories, and Bean and Thayer | 2Jof the National 
Bureau of Standards, establishes that the lapse of refrac- 
tive index, dn/dz, where n is the refractive index and z 
the heieht. above. the earth, .is not. constant, as it was 
assumed to =e in the 4/3-earth model. These investira- 
Goms suceest that. the, revised standard n=profile. i.s.more 
suitably described in the mean by an exponential curve, 
Father than i9,-linear one,.,.such that: 


~G# (1) 


ne il+ce 
where €; and Co are empirically determined constants, 

At this point, let us consider, from a quantitative 
viewpoint, the effiect-of the index of refraction on a 
radar ray. Since n is generally of the order of 1.0003, 
it is more convenient to discuss its scaled-up counter- 
part. refractivity, N, wheres 

N=(n-1)10° (2) 
Phereftore, if n=1.005, the refractivity bs N=300, En 
the mean, the vertical eradient? of N for the lower layers 
of the atmosphere is 12 units/1l000 feet. Substituting 


this value into the relative curvature eovation [ see Ref. 2] 


Ifhe vertical gradient will be represented by — 7 > 
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with the average or ' 


'standard" curvature, resulting from 
the 4/Z3-earth model. For -S2 = 47.8 units/1000 feet, 
Equation (2)shows that the radar-ray curvature is the same 
a5s.that.of.the. earth. This, condition, is .known as. a- trap= 
Dung, or ducting. condition. Sfnce vertical Negradients 
well in excess of 48 units/1000% feet exist in the atmos- 
phere, the radar-ray path.can be.significantly altered 
during its passage through the atmosphere. .To illustrate 


bae error. this. bending:.introduces.in altitude determining 


radars, consider a.transmitter,T situated below a super- 


perractive layer .in Figure 2. 


The solid line TS, is the actual ray path between the 
target and the transmitter. The dashed line, TS', is the 
ray path along which one assumes the radar energy has travel- 


ed, if the 4/3-earth model applies, and from which the s 


oO 
cr 


computes its displayed height. Therefore, the radar set 
will display a height H', which, as can be seen, constitutes 


a Sizeablée departure from the true height, H. 
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An exponential lapse of n with heizht has been assumed 
both below and above the trapping layer, followin Sovation 
(1). A single trapping layer has been assumed, of thickness 
Aiotmnd exhibiting athinearuilapse of hn.s As can be iseen 
from the eeometry of Figure 5 and Equation (1), the value 
of n at any point below the base, Z=h, of a trapping layer 
Zesexpress.ed ¢ 


—C,2 
ee ak 


Similarly, the value of n at any point within the trapving 
layer of thickness Ah can be expressed as: 
een ANY E 
oo eee )+ ay (z-h) hgzehrah (5) 
Finally, the value of n above the trapping layer can be 


expressed; 


cemarge Ce 2 ou en Gh [ereah ey af an ah 
Note that the model just constructed has the same slope 
dn/dz above the trapping layer as it would have had in 
the absence of a trapping layer. Whether this is fulfilled 
in practice is not known; however, this feature of the 
model is very convenient from a computational viewpoint. 

This feature may be stated another way: for any given 
trapping layer of specified thickness, gradient and height, 
Equation @) can be seen to have the form 

ne 4+Ce%* a Ki) (7) 

where K(t) is a function of the properties of the trap- 


ing layer, and is a constant for any set of computations 
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dealing with any single trapping layer. 
Faving described the model, the next step is the 
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selection of appropriate values for Cz and Coe The 
gations of Bean and Thayer | 2] have established a value for 
the constants based on 39 station-yesns of radiosonde 
observations, representing an all-season average atids 
climatically diverse locations. On the other hand, Bawer, 
Mason and Wilson Li] have Sugzested values for the constants 
based on a selected set of soundings representing a weooL, 
dry" spring atmosphere. While the work of Bean and Thayer 
represents a more precise representation of an average "all 
cases" sounding, the atuosphere described by Bauer, kkason, 
and Wilson better approximates an ideal atmosphere into 
which anomolous conditions can be interjected. Therefore, 
the values proposed by Bauer, Mason, and Wilson have been 
used to define the exponential atmosphere, and Equation (1) 


now becomes; p 
=<O3 TIS 
m=. 1+ ,00032.e ( 


(e9) 


where 2 is measured in kilofeet. 
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Figure 4. Ray geometry. 

Consider the ray in Figure’4 entering a spherically 
homogeneous layer of the atmosphere with entrance angle 6; ; 
and emerging from the top of the layer with anple of emer- 
gence Ay, ; Angeles By and Bis; are always considered to be 
measured from the local horizontal. Snell's law tells us 
[3]: 

hi; (a+z;)cos B ee, (a + Zu) cos Aévt = Constant (9) 
and nj,q ars the indices of refraction at tottom 
and top of this layer respectively; 2; andZj;,, are the 
heights of the points above the surface of the earth; 


3 2 n 2 (S) ° 
and a is the true radius of the earth*. Expressing (9) 


2. J s 2 
“For the purposes/of this paper, ythe earth is assumed a 
perfect sphere of radius 20898.908 kilofeet. 
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the last equation may be written: 
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Replacing the differentials with a finite-difference fom, 


(10) can be written: 


AS: = Bos a Bi Selo. 
Cos B [ —2Chia-ni) Es 1 (11) 
(Riss #1; Gia 72) A+ (Zt%us 


| s Rr 7} 5 GC . 5 og 
Using Eouation (¢) to solve for ,,,-8;, one obtains: 


ee ia = hi (a+ 2%) CosB 5 
fis Pr - ( Mira Ot Bea) A) ~ pi aa 


Substituting (12) inte (11) vives: 


feat ni (arzi) Cos Bi | - B, 


ied ( Qt Eur: a) 


AS; = 
| ral Nisa 7 A } i 1 
Cos Bi d u) 


(13) 


(Neg tM Meme) Ot (2:42.44) 
2 


Houation (13) can be numerically solved by the following 
steps: 
1. select small increments AZ; =2,,,-2; , given an 
initial angle fi 
Za compute AS; and Bisa over small thicknesses using 
Houations (4), (5), (6),(12). and (13) 


S. compute range S and height 4 where: 


o= > Asi - = S Ever (14) 
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For computations used in this paper, AZ; 
200 Peet, 

Let us look for a moment at Equation (12). If nj 
is aufsvetomt ereater than ny, , and the ray has a low 
angle of elevation so that Cosni; 24 , the tem Be Gos Bi 
Will«abecomesagreater.thanounity, ».Under+these.conditions the 
Gopune (of an,angle «cannot exist; .and.Snell's;law’ failss 
smehaancondi tion. represents«ceomplete reflection of »therray, 


Srrerepping.,.Under,»these «conditions, -+the-heitght.of pany 


bars 


eun.cetbected by, Chis: ray. widl,bepless «than, the,heisht 


Ghrthentop.of,the trapping layer. 
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4, Computation procedure 

The model as described in Sections. “was computed fer 
Ve top of “the ‘trapping Layer "at sévery "L000 feet up to 
5000 feet, and for thicknesses of 200, 400 and 800 feet, 
Wwiuy N-gradients of 50, 100, and 200 units/1000 feet. 
These values have been selected so that all cormonly- 
ecturring atmospheric trapping layers will lie within the 
ranpe Sot /~paraméter-valués “selected, “providing -a-basis “for 
pocular interpolation.” Sélection of thé numberof parameters 
used has béen a compromise between compléte tabular listings 
Bad bine SiPr7é-of Appendix 2.8 “The ‘author Has “tried to select 
be cases so that both requirements are adequately met, 
keeping Appendix 2 to such a size that its value will not 
Bemiést.in its bulk. 

The procedure outlined in Section 5S was programmed 
for the CDC-16C4 high speed digital computer. Rays were 
computed for increments in the initial angle fo of one- 
eighth degree up to one-half degree; then in half-degree 
increments from 6 to 20 degrees, then 5 degree increments 
for angles from 25 to 60 degrees. Heights were computed 
in200 feet steps, and for each height the range computed. 
The computed range and the known initial angles were then 


used to solve the 4/3 earth heicht using the Equation (1): 


a cr 
z=S ee B. + 352 (15) 


where @ is the effective earth's radius, in this case (4/3)a. 
Equation (15) is noted to -be an approximation formule, but- 
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it .zonstitutes the.commonly employed approximation formule 
in the design of .altitude-determining radars., Since,this 
paper concerns corrections to radar observaticns, rather 
than corrections to the 4/3-earth model, the approximation 
is appropriate. 

It will be noted that for each simultaneous computation 
of a complete ray through the model atmosphere and the 4/3 
earth atmosphere, both rays have a common initial angle and 
range. The common range value for the two rays was reduced 
PoOmioGLenvals of «ten miles,..and the corresponding reishts 
for each model were interpolated to proper values for the 
10-mile intervals for identical values of B, . After all 
rays have been computed, two sets of tables will exist for 
@ach range, one set containing a family of heishts for the 
4/3-earth atmosphere, and the second table containing corre- 
sponding values of the model height. The table containing 
the 4/3-earth model heipht was then reduced to 10,000 feet 
intervals, and the corresponding true heights were linearly 
interpolated to their appropriate values. The general 
technique can be described in general terms as two-way 
tabular interpolation, employing the search and logic functions 
ef the CDC-1¢€04 computer’ The results were then printed in 
the form displayed in Tables 1-45 


For cases where ray trapping occurred, the computer 


Ome author has maintained a complete set of the IBN cards 
used in the main program, ame veple to those interested in 
the precise programming techniques. 
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generated ranges to the point of trapping and these were 


placed in a range table opposite the 4/3-esarth heicnt. 


Heights were obtained with 
model heiechts. The actual 
arbitrarily defined as the 
surrace’ and the -top ‘ot the 
computations were made for 


after it reached the level 


their corresponding 4/3-earth 
height of any trapped ray was 
mid-point height between the 
Grapping layer.’ "No further 

any ray in the model atmosphere 


of trapping Vas just ‘delined:. 


mre ray corresponding toy the next Increment’ of fo was then 


computed as outlined in Sections 3 and 4. 
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De Conclusions. / 

Test runs of the computation scheme for a strict 
exponential model were in close agreement with the publish- 
ed results of Bauer, Mason, and Wilson, being. within? 20C 
Peet» Ofte these rresults,, Whterttables: in- Appendix P-should 
therefore, provide. adequate accuracy for most radar uses. 

At» this point, the author would like to caution the 
mepoer that, while the digital computer provided results to 
eight cecimal places, at extreme ranges and low altitudes, 
the tables probably do not provide accuracy to any greater 


——— 


degree thant 500 feet. In order to provide greater accuracy, 


smaller height increments and more sophisticated integration 
procedures would have to be used. 

The tabulated values in Appendix I should provide a 
valuable addition to the height-finding radars. The ex- 
istence of trapping layers can be ascertained using pro- 
cedures outlined in various AROWA publications U5]. 

This problem constituted the first computational pro- 
gram of any sige completed on any of the CDC-1604 computers. 
As a result, the task was undertaken prior to the existence 
of an adeouate subroutine library for the computer. Since 
the author not only had to program the primary routine, as 
well as all subroutines, an approach was selected to keep 
these subroutines to a minimum. For example, a natural loga- 
rithm subroutine and other types of interpolation subroutines 
could have been useful, but limited time precluded the develop- 


ment of such routines during this particular study. 
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APPENDIX 1 
Use of Tabies 

Upon selecting the proper table for use, based on pre- 
vious knowledge of atmospheric conditions enter the left 
hand column with the observed range, and enter across the 
top with the observed radar height, in thousands of feet. 
The true altitude of the target may then be read from the 
table. 

EXAMPLE: Assume we are using Table 1, with an ob- 

served height of 20,000 feet (20 kilofeet) and an 

observed ranzecof 150 mites: ;s the wsable gives: 

+,20612005 +02, Note,that thernvalue is listed as a 

fractional coefficient with any accompanying power 

of ten. Moving the dectmea point two places to the 
right as indicated by the +02, the value becomes 

20612005 kilefeet, or 20,612 feet. Considering tie 

accuracy of the computations, this would te recorded 

as) £0,500 fee, 

In the event trapping eecunred. (a notation! has been 
ineéluded at. the bottom. of the table indicating the first 
angle, of those angles tested, which passed through the 
layer without being trapped. Elevation angles less than 
this angle are to be treated with caution, as the height 


values in the table will indicate. 


luxcept if trapping does not eecur for any initial ray 
BA) 


angle; then no notation appears under the table. 
15 
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APPENDIX 2 


TABULAR RANCE-HEIGHT VALUES 


Tables 1-45 
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HEIGHT DETERMINATION FOR MODEL ATMOSPHERE WITH TRAPPING LAYERS 
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HETGHT DETERMINATION FOR MODEL ATMOSPHERE WITH TRAPPING LAYERS 
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HEIGHT DETERMINATION FOR MODEL ATMOSPHERE WITH TRAPPING LAYERS 
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HEIGHT DETERMINATION FOR MODEL ATMOSPHERE WITH TRAPPING LAYERS 
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